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Abstract 



We investigate the fermionic condensate and the vacuum expectation value of the energy- 



momentum tensor for a massive fermionic field in the geometry of two parallel plate on the 
background of Minkowski spacetime with an arbitrary number of toroidally compactihed 
£vq \ spatial dimensions, in the presence of a constant gauge field. Bag boundary conditions are 

■ imposed on the plates and periodicity conditions with arbitrary phases are considered along 

the compact dimensions. The nontrivial topology of the background spacetime leads to 
an Aharonov-Bohm effect for the vacuum expectation values induced by the gauge field. 
The fermionic condensate and the expectation value of the energy-momentum tensor are 
periodic functions of the magnetic flux with period equal to the flux quantum. The boundary 
induced parts in the fermionic condensate and the vacuum energy density are negative, 
with independence of the phases in the periodicity conditions and of the value of the gauge 
potential. Interaction forces between the plates are thus always attractive. 

However, in physical situations where the quantum field is confined to the region between 
the plates, the pure topological part contributes as well, and then the resulting force can be 
either attractive or repulsive, depending on the specific phases encoded in the periodicity 
conditions along the compact dimensions, and on the gauge potential, too. Applications of 
the general formulas to cylindrical carbon nanotubes are considered, within the framework of 
a Dirac-like theory for the electronic states in graphene. In the absence of a magnetic flux, the 
energy density for semiconducting nanotubes is always negative. For metallic nanotubes the 
energy density is positive for long tubes and negative for short ones. The resulting Casimir 
forces acting on the edges of the nanotube are attractive for short tubes with independence 
of the tube chirality. The sign of the force for long nanotubes can be controlled by tuning the 
magnetic flux. This opens the way to the design of efficient actuators driven by the Casimir 
force at the nanoscale. 



*Also at Tomsk State Pedagogical University, Tomsk 



1 



1 Introduction 



In a good number of problems one needs to consider a physical model on the background of 
some manifold with compactified spatial dimensions. Many high-energy theories of fundamental 
physics are formulated in a higher dimensional spacetime and it is commonly assumed that 
the extra dimensions are compactified. In particular, additional compact dimensions have been 
extensively used in supergravity and superstring theories. Prom the inflationary point of view, 
universes with compact dimensions, under certain conditions, should be considered a rule rather 
than an exception [lj. Models of a compact universe with non-trivial topology may play an 
important role by providing proper initial conditions for inflation. There has been a large activity 
to search for signatures of non-trivial topology by identifying ghost images of galaxies, clusters 
or quasars. Recent progress in observations of the cosmic microwave background provides an 
alternative way to observe the topology of the universe. An interesting application of the field 
theoretical models with compact dimensions recently appeared in nanophysics [2]. The long- 
wavelength description of the electronic states in graphene can be formulated in terms of the 
Dirac-like theory in 3-dimensional spacetime with the Fermi velocity playing the role of speed 
of light (see, e.g., Refs. OHj). Single-walled carbon nanotubes are generated by rolling up a 
graphene sheet to form a cylinder and the background spacetime for the corresponding Dirac-like 
theory has topology R 2 x S . For another class of graphene-made structures, called toroidal 
carbon nanotubes, one has as background topology R 1 x (5 1 ) 2 . 

The boundary conditions imposed on fields along compact dimensions give rise to a modifica- 
tion of the spectrum of the vacuum fluctuations and, as a result, to Casimir-type contributions 
in the vacuum expectation values of physical observables (for the topological Casimir effect 
and its role in cosmology see [5 J and references therein). In models of Kaluza-Klein type, the 
Casimir effect has been used as a stabilization mechanism for moduli fields and as a source for 
dynamical compactification of the extra dimensions, in particular, for quantum Kaluza-Klein 
gravity (see Ref. [6]). The Casimir energy can also serve as a model for dark energy needed 
for the explanation of the present accelerated expansion of the universe (see [7] and references 
therein). In addition, recent measurements of the Casimir forces between macroscopic bodies 
provide a sensitive test for constraining the parameters of long-range interactions, as predicted 
by modern unification theories of fundamental interactions [8j. The influence of extra compact- 
ified dimensions on the Casimir effect in the classical configuration of two parallel plates has 
been recently discussed in [9J, for the case of a scalar field, and in [ID], for the electromagnetic 
field with perfectly conducting boundary conditions. 

More recently, interest has been focussed on the topic of the Casimir effect in braneworld 
models with large extra dimensions. This type of models (for a review see |llj ) naturally appear 
in the string/M theory context and they provide a novel set up for discussing phenomenological 
and cosmological issues related with extra dimensions. In braneworld models the investigation 
of quantum effects is of considerable phenomenological value, both in particle physics and in 
cosmology. The braneworld corresponds to a manifold with boundaries and the bulk fields 
will give Casimir-type contributions to the vacuum energy and, as a result, to the vacuum 
forces acting on the branes. The Casimir forces provide a natural mechanism for stabilizing 
the radion field in the Randall-Sundrum model, as required for a complete solution of the 
hierarchy problem. In addition, the Casimir energy gives a contribution to both the brane and 
the bulk cosmological constants. Hence, it has to be taken into account in any self-consistent 
formulation of the braneworld dynamics. The Casimir energy and corresponding Casimir forces 
within the framework of the Randall-Sundrum braneworld [12] have been evaluated in Refs. [13] 
by using both dimensional and zeta function regularization methods. Local Casimir densities 
were considered in Ref. p3]. The Casimir effect in higher dimensional generalizations of the 
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Randall- Sundrum model with compact internal spaces has been investigated in |15j . 

In the present paper we will study the fermionic condensate, the Casimir energy density 
and the vacuum stresses for a massive fermion field in the geometry of two parallel plates on 
a spacetime with an arbitrary number of toroidally compactified spatial dimensions. We will 
impose generalized periodicity conditions along the compact dimensions with arbitrary phases 
and MIT bag boundary conditions on the plates. The presence of a constant gauge field will 
be assumed as well. Though the corresponding field strength vanishes, the nontrivial topology 
of the background spacetime leads to Aharonov-Bohm-like effects on the vacuum expectation 
values. The total Casimir energy in the geometry under consideration has been discussed in 
Ref. [16], in the absence of a gauge field. The investigation of local physical characteristics in the 
Casimir effect, such as the expectation value of the energy- momentum tensor and the fermionic 
condensate, is of considerable interest. Indeed, local quantities contain more information on the 
vacuum fluctuations than global ones. In addition to describing the physical structure of the 
quantum field at a given point, the energy-momentum tensor acts as the source in Einstein's 
equations and, therefore, it plays an important role in modeling a self-consistent dynamics which 
involves the gravitational field. The fermionic condensate plays an important role in the models 
of dynamical chiral symmetry breaking (see review [17] for chiral symmetry breaking in the 
Nambu-Jona-Lasino and Gross-Neveu models on the background of a curved spacetime with 
non-trivial topology and |18j for very recent developments). 

The fermion Casimir energy for two parallel plates in 4-dimensional Minkowski spacetime 
with trivial topology has been considered in [19], for a massless field, and in |20| in the massive 
case. For arbitrary number of dimensions, the corresponding results are generalized in Refs. |21[ 
[22] for the massless and massive cases, respectively. The fermionic condensate for a massless 
field has been considered in Refs. |23j . The Casimir problem for fermions coupled to a static 
background field in one spatial dimension is investigated in [23]. The interaction energy density 
and the corresponding force are computed in the limit that the background becomes concentrated 
at two points. The fermionic Casimir effect for parallel plates with imperfect bag boundary 
conditions modelled by J-like potentials is studied in [25]. The topological Casimir effect and 
the vacuum expectation value of the fermionic current for a massive fermionic field in a spacetime 
with an arbitrary number of toroidally compactified spatial dimensions have been considered in 

[26J W\. 

The paper is organized as follows. In the next section, we specify the eigenfunctions and 
the eigenmodes for the Dirac equation in the region between the plates assuming bag boundary 
conditions on them. The fermionic condensate in this region is considered in Sect. [3l By 
using an Abel-Plana-type summation formula, we express the condensate as the sum of a pure 
topological, single plate contribution and interference parts. Various limiting cases are then 
considered. The vacuum expectation value of the energy-momentum tensor is investigated in 
Sect. HI A generalization for a conformally-flat background spacetime is given. In Sect. [5] 
we study applications of the general formulas to the Casimir effect for electrons in a carbon 
nanotube, within the framework of a 3-dimensional Dirac-like model. The main results of the 
paper are summarized in Sect. [6] 

2 Fermionic eigenfunctions 

Consider a spinor field, if), propagating on a (D + l)-dimensional flat spacetime with spatial 
topology R p+1 x (5* 1 ) 9 , p + q + 1 = D. We will denote by z p+ i = (zi, . . . ,z p+ \ = z) and 
z q = ( z p+2i ■ ■ ■ i z d) the Cartesian coordinates along the uncompactified and the compactified 
dimensions, respectively. For these coordinates we have — oo < z\ < oo, I = 1, . . . ,p + 1, and 
^ z% ^ Li for I = p + 2, . . . , D, with Li being the length of the l-th compact dimension. We 
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assume that along the compactified dimensions the field obeys quasiperiodic boundary conditions 

rP{t,z p+1 ,z q + Liei) = e 27Tiai ^(t,z p+1 ,z q ), (2.1) 

with constant phases |orj| 1/2 and e; is the unit vector along the direction of the coordinate 
Zi, I = p + 2, . . . , D. Condition (|2.1|) includes the periodicity conditions for both untwisted and 
twisted fermionic fields as special cases with a\ = and a\ = 1/2, respectively. The special 
cases ai = 0, ±1/3 are realized in nanotubes. 

In this paper we are interested in the fermionic condensate and in the vacuum expectation 
value (VEV) of the energy-momentum tensor induced by two parallel plates located at z = 
and z = a. On the boundaries the field obeys MIT bag boundary conditions 

(l + iyn M )V = 0, 2 = 0, a, (2.2) 

with 7^ being the Dirac matrices and the outward oriented (with respect to the region under 
consideration) normal to the boundary. Note that from the conditions (|2.2p it follows that, on 
the boundaries, vptp = and n^j^ip = 0, where ^ = ifj^ is the Dirac adjoint and the dagger 
denotes Hermitian conjugation. In the discussion below the calculations will be done for the 
region between the plates, < z < a, where we have = — 5fi +1 at z = and = <5 P +1 at 
z = a. The expressions for the VEVs in the regions z < and z > a are obtained as limiting 
cases. 

Dynamics of the massive spinor field is governed by the Dirac equation 

vfdpil) - mi/) = . (2.3) 

In the (D + l)-dimensional spacetime, the Dirac matrices are Nd x No matrices with Nd = 
2 [(£>+i)/2] ) where the 

square brackets mean integer part of the enclosed expression. We will take 
these matrices in the Dirac representation: 

From the anticommutation relations for the Dirac matrices one has cr^a^ +ov<7+ = 25 In the 
case D = 2 we have Nd = 2 and the Dirac matrices are 7^ = (crp3, icrpi, icp2), with crp^ being 
the 2x2 Pauli matrices. 

The boundary conditions (|2.ip and (|2.2jl lead to the modification of the spectrum for vac- 
uum fluctuations of the fermionic field and, as a result, to the topological and boundary induced 
Casimir effects on the VEVs of physical observables. For the evaluation of the VEVs we need 
the complete set of positive- and negative-energy solutions to the Dirac equation satisfying the 
boundary conditions (12. 2p . The dependence on the coordinates parallel to the plates, Z11 = 
(21, ... , z p , z p+ 2, • • • , zd), can be presented in the standard exponential form exp(«kii • zm), with 
k|| = (kp,k 9 ) and k p = (k%, . . . , k p ), k ? = (k p+ 2, ■ ■ ■ , &d). The eigenvalues for the compo- 
nents of the wave vector along the compactified dimensions are determined from the periodicity 
conditions (|2 . 1 [) : 

k q = (27r(n p+2 + a p+2 )/L p+2 , . . . ,2n(n D + a D )/L D ), (2.5) 

with n p +2> ■ ■ ■ j no = 0, ±1, ±2, .... For the components along the uncompactified dimensions one 
has —00 < ki < 00, I = 1, . . . ,p. The corresponding positive- and negative-energy eigenspinors 
have the form 

*P = A^( ia - Vx, ^ + m) Y (2.6) 
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where cr = . . . , <td), qj = yk| + + k^ + m 2 and /? is the collective index for the set of 

quantum numbers specifying the solutions (see below). The spinors in (|2.6p are given by the 
expressions 

ik ir z ll (in , P ik P+i- z + in P ~ ik P+i z 



ip = e Jk ir z ll lip + e lk p+ iz + tp_e 



X = e" ik l|- z l 



(x+e ikp+lZ + x-e~ ikp+iz ) , (2.7) 



From the boundary condition (|2.2I) on the plate at z = we find the following relations 
between the spinors in (|2.7p 

m(oj + m) + - /c p+ i<7 p+1 crt • ky 
^+ = ? ^ w — , \ V-> 

m(uj + m) + £^ +1 - /c p+1 ct+ +1< T|| • ky 

X- = 7 r-r, 77 — : \ X+, (2.8) 

(m + «fc p+ i)(u; + mj 

where crii = (a\, . . . , a p , a p +2, ■ ■ ■ , &d)- We will assume that they are normalized in accordance 
with ipt-tp- = X+X+ = 1- As a set of independent spinors we will take <£_ = vM' and x+ = w^' , 
where w/ CT \ a = 1, . . . ,Nd/2, are one-column matrices having N^/2 rows with the elements 



w 



and u/ CT )' = iw^ a \ Now the set of quantum numbers specifying the eigenfunctions 
(j2.6[) is /3 = (k, cr). From the boundary condition at z p+1 = a it follows that the eigenvalues of 
k p+ i are roots of the transcendental equation 

masm(k p+ ia) / (k p+ ia) + cos(/c p+ ia) = 0. (2-9) 

All these roots are real. We will denote the positive solutions of Eq. (|2.9[) by X n = k p+ \a, 
n = 1,2,.... For a massless field one has A n = 7r(n — 1/2). Note that Eq. (I2.9h does not contain 
the parameters of the compact subspace and is the same as in the corresponding problem on 
the topologically trivial Minkowski spacetime (see [5]). Note that this will not be the case in a 
more general class of compact subspaces. 

The normalization coefficient Ap in (|2.6p is determined from the orthonormalization condi- 
tion 

dz|| £ dz p+1 ^ ±)+ ^t } = 5pp. (2.10) 

Here, the symbol 5pp> is understood as the Dirac delta function for continuous indices and the 
Kronecker delta for discrete ones. The substitution of the eigenfunctions (j2.6|) into this condition 
leads to the result 

sm{2k p+ ia 



2 oj + m 
13 = 4(27r)PwaV r 5 



1 



2fc p _|_ia 



(2.11) 



where V q = L p+ 2 • • • Ljj is the volume of the compact subspace. 

We can generalize the eigenfunctions given above to the situation when an external elec- 
tromagnetic field with vector potential A p = const is present. In spite of the fact that the 
corresponding magnetic field strength vanishes, the non-trivial topology of the background 
spacetime leads to the appearance of an Aharonov-Bohm-like effect for the physical observ- 
ables. In particular, the corresponding VEVs depend on A^. Now the Dirac equation has the 
form ^^(d^+ieA^ip — mip = and, by making use the gauge transformation A^ = J 4^+9 At A(x), 
tp(x) = ifj' {x)e~ te ^ x \ with the function A(x) = A^x^, we see that the new function il>'{x) sat- 
isfies the Dirac equation with A'^ = and the quasiperiodicity conditions similar to (|2.ip with 
the replacement 

ai -> ai = ai + eA l L l /{2i:). (2.12) 
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The eigenvalues for the wave vector components along compact dimensions are defined by ki = 
2ir(ni + ai)/L[ and the corresponding eigenspinors are obtained from those given above with the 
replacement (|2.12|) . 

3 Fermionic condensate 

The fermionic condensate is among the most important quantities that characterize the prop- 
erties of the quantum vacuum. Although the corresponding operator is local, due to the global 
nature of the vacuum, this quantity carries important information about the global properties 
of the background spacetime. Having the complete set of eigenspinors, we can evaluate the 
fermionic condensate by using the mode-sum 

(^} = ^4- ) (x)^ ) (x), (3.1) 

P 

where (• • • } stands for VEV. By taking into account the expression (12. 6p for the negative-energy 
eigenspinors, the fermionic condensate in the region between the plates can be expressed as 

Nd I rfkp sin(A n z/a) 



q n q eZi J y ' n=\ 

^ ma sm(X n z/a) + X n cos(X n z/a) 

, 1-J.2J 



1 -sin(2A n )/(2A 
where n q = (n p+ 2, • • • , n£>) and 



u = / a 2 + k 2 + k 2 + m 2 . (3.3) 

Of course, the expression on the rhs of Eq. (|3.2p is divergent. We will assume that some cutoff 
function is present, without writing it explicitly. 

For the further evaluation of the fermionic condensate we apply to the sum over n in Eq. (|3.2p 
the Abel-Plana-type summation formula 

^ *HK) nmam [°° . [°° f(ix) - fj-ix) 

g i-s m( 2A„)/ ( 2A„r -2(^TT) + i */w-y * ' (3 - 4) 



n=l 

with the function 



f(x) = — ; = = = = [ma sm(x z / a) + xcos(xz/ a)]. (3-5) 
c 2 + k 2 a 2 + ml a 2 



Formula (|3.4p is a special case of the summation formula derived in Ref. [29] on the basis of the 
generalized Abel-Plana formula (see also Ref. [30]). In Eq. (13. 5p and the discussion below we 
use the notation 

D 

ml q = k 2 nq + m\ k 2 nq = [M"l + <*l)/Li} 2 . (3.6) 

l=p+2 

After the application of formula (|3.4p . the fermionic condensate is split into 

= W} (0) + W) (1) -^ E fw^Tr dx 



sinh(xz) msinh(xz) + xcosh(xz) 



9l o 9 x+m 2ax _i_ i 

"r — Kt — mz. x-m c ~ 

V n q 



(3.7) 
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n 9 eZ^ ^ Jk2 +1+m 2 ? 



where 



is the fermionic condensate in the topology R p+1 x (S" 1 ) 9 when the boundaries are absent. The 
term 

Njj ^-^ f dkp f°° x s'm(2xz) — mcos(2xz) 



is the part induced by the plate at z = when the second plate is absent. The last term on the 
right of formula (j3.T|) comes from the last term in Eq. (|3.4h and it is induced by the presence of 
the second plate. Note that this term vanishes at z = 0. 

For points away from the boundaries, the boundary induced part is finite and the cutoff 
function in the corresponding expressions can be safely removed. Renormalization is needed for 
the purely topological part only. The latter has been investigated in Ref. [26] (for the topological 
fermionic Casimir effect in de Sitter spacetime with toroidally compactified spatial dimensions 
see Ref. [31] ). By making use of the zeta function technique, the corresponding renormalized 
expression is presented in the form 



W) (0) = - ( 27T )( D D+m cos ( 27rn g ' a i 

V ' n,€Z9 



f(D-i)/2(™^L 2 p+2 n 2 p+2 + ■■■ + L 2 D n 2 D ) 
(Ll +2 n 2 p+2 + ---+L 2 D n 2 D )^)/ 2 



(3.10) 



with a q = (a p+ 2, ■ ■ ■ , «d), and 



/„(*) = z u K v [z). (3.11) 



The prime on the summation sign in (|3.10p means that the term n q = is excluded from the sum. 
An alternative expression for {^pip)^ is obtained in Ref. [26] using the Abel-Plana summation 
formula. In the discussion below we will be concentrated on the boundary induced parts. 

For further transformation of the single plate part (j3.9l) , we write the function in the integrand 

as 

x sin(2xz) - m cos(2x2:) = — - [(m + ix) e 2lxz + (m - ix) e~ 2txz ] . (3.12) 

In the integral over x in (|3.9p we rotate the integration contour by an angle 7r/2, for the term 
with the exponent e 2lxz , and by — 7r/2, for the term with the exponent e~ 2lxz . As a result, we 
get 

It follows from this expression that (^ip)^ is always negative. By using the relation 



f(x)dx 7t(p+i)/2 r°° 

. J v ' = = — r—- I dx(x 2 -ml ) {p ~ 1)/2 f(x), (3.14) 

^f^< q Jx 2 -k 2 - ml q T ((P + l )l 2 ) Jmn t 



we find 



^)W = MEJ2 r dx(x 2 -m 2 nq )^/ 2 (m-x)e- 2xz , (3.15) 
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A P = (3-16) 



with the notation 

(47r) -( P +l)/2 

r((p + i)/2) 

The integral in Eq. (|3.15p is expressed in terms of the modified Bessel function of the second 
type, namely 



= N { ^y/2+iV q E [2mzf p/2 (2m nq z) - f p/2+1 (2m nq z)] . (3.17) 



In the absence of compact dimensions, from Eq. (I3.17P one finds 

W>S> = (^) ( (D+iV2 [2m*/(D-i)/ 2 (2mz) - / (D+ i )/2 (2w)] . (3.18) 

Let us consider some limiting cases of the general formula (|3.17p . In the limit when the 
length of one of the compactified dimensions, say z J , j ^ p + 2, is large, Lj — > oo, the dominant 
contribution to the sum over rij in Eq. (|3.17|) comes from large values of rij and we can replace 
the corresponding summ by an integral, with the help of the relation 

+oo 

7T 



poo 

V f{2Tr\n j + a j \/L j )^ / dy/fo). (3.19) 
The integral is evaluated by using the formula 



J n 3 - 



1 pK 



dy fv{c\Jy 2 + b 2 ) = - J -fv+i/2(bc), 

and we can see that, from Eq. (|3.17p . the corresponding formula is obtained for the topology 
R p+2 x (S" 1 ) 9-1 . At small distances from the boundary, z <C m _1 ,L;, the main contribution to 
the series in Eq. (|3.17p comes from large values of ni and we can replace the summation by the 
integration and to the leading order we find 

This leading behavior does not depend on the lengths of the compact dimensions and, as it is 
seen from Eq. (I3.18p . coincides with boundary induced part of the fermionic condensate for a 
single plate in a space with trivial topology R D in the case of a massless field. 

Now, let us consider the limit Li <C z. In this case, and for ai = 0, the main contribution 
comes from the zero mode with n q = and to the leading order we find 



mz) 



Comparing with ([3. 18ft . we see that the quantity Vq^ip)^ / N £, coincides with the corresponding 
result for a plate in topologically trivial (p + l)-dimensional space, R p+1 . The contribution of the 
nonzero modes is exponentially suppressed and, for ai ^ 0, the zero mode is absent. Assuming 
that mz is fixed, to leading order we have 

(frANfi) - ^D< +1 e~ 2moZ (321) 
m) - 2V q (4nm z)(P+W ^ } 
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where 

D 

ml = (^m/Li) 2 . (3.22) 

l=p+2 

In this case, the boundary induced part in the fermionic condensate is exponentially suppressed. 

Using (|3.14p . we can also simplify the expression for the second plate induced part in 
Eq. (|3.7p . Combining with Eq. A3. 15f> we find 

,t a /7/ \(0) A p N D ^ r 4 (x 2 -< q ) {p - 1)/2 

x [(m + x)(e 2xz + e 2ax - 2xz ) - 2m] , (3.23) 

where the second term on the rhs is the boundary induced part, which is always negative. For 
a massless field, by using the expansion (e y + 1) _1 = — X)nLi( — l) n e~ ny , from this formula we 
find 

V ' q n„eZi n=l 



\- K p/2+1 (2k nq (an - \ aj - z\)) 



where a% = and 02 = a and where k nq is defined by Eq. (|3.6p . In the case of trivial topology 
R D , from here we obtain 



- N D T((D + l)/2) ™ 

v ' n=l 



+ 



{n — z/a) D (n—l + z/a) 



D 



(3.25) 



The series in these formula are given in terms of the Hurwitz zeta function. The last expression 
can be further simplified, for odd values of D, as 

- _ N D T((D + l)/2) d - 1 1 

{W)BP - - 2D+ 2 7T (D-l)/2 a D dx D-l s[n{7Tx) x=z/a ( 3 - 26 ) 

In the special case D = 3, from here we obtain the result given in Ref. |23j . 

Extracting the parts corresponding to the single plates, the fermionic condensate can also 
be presented in the form 

= (^) (0) + E (^)f + A ^)' ( 3 - 27 ) 

where ('tp^j is that part in the condensate induced by a single plate located at z = dj, with 
a\ = and 02 = a, and the interference term is given by the expression 

A P N D ^ p , (x 2 -<) (p - 1)/2 



A(^) = ^ £ r dx 



x+m, c 2aa; _|_ j 



x—m 



x [2m + (s - m)(e- 2a;z + e 2 ^" 2 ^)] . (3.28) 

As it is seen from this formula, the interference part in the fermionic condensate is always 
positive. Note that the divergences on the plates are contained in the single plate parts and that 
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the interference term is finite for all values ^ z ^ a. For a massless field, similar to Eq. (|3.24[) . 
one finds the expression 

A <**> = -™f E < 2+i D-ir E <**» 

v ; g n q eZi n=l 3=1,2 v 1 ■> w 

In the limit Lj<a and for az = 0, the main contribution comes from the zero mode with = 
and to leading order we find A(tpip) « N£)A(tptp) p> P +i / (N p+ \V q ), where 



CO 



( x 2 _ m 2)(p-l)/2 



A(^)bp+i = A^P+l / da; B+m „ — — 

x [2m + (s - m)(e" 2 ^ + e 2xz - 2ax )] , (3.30) 

is the interference part in the fermionic condensate for the geometry of two parallel plates in 
(p + l)-dimensional spacetime with trivial topology RP +1 . When a\ ^ 0, there is no zero mode 
and the interference part is exponentially suppressed by the factor e~ 2am ° , where mo is defined 
in Eq. 4222]). 

In the discussion above we have considered the fermionic condensate in the region between 
the plates, < z < a. For the regions z < and z > a, the expression for the fermionic 
condensate has the form 

(^/>) = (^)(°) + (3.31) 

where the plate-induced part is given by Eq. (|3.17p with the replacement z — > \z\, in the region 
z < 0, and with the replacement z — > z — a, in the region z > a. 

In the presence of a constant gauge field A^, the corresponding formulas for the fermionic 
condensate are obtained from those given above by the replacement (|2,12p . From these formulas 
it follows that the fermionic condensate is a periodic function of A[Li with the period of the flux 
quantum <I>o = 2vr/|e| (^o = 2Trhc/\e\ in standard units). It is an even function of a/. 

The formulas corresponding to the special case with a single compact dimension are obtained 
from those given above by taking p = D — 2, q = 1, and k Uq = 2-K\riD + od|/-Ld. In Fig. [T]we 
plot, for the simplest Kaluza-Klein-type model with spatial topology R 3 x 5 1 , the dependence 
of the fermionic condensate for a massless fermionic field vs the ratio z/a, in the region between 
the plates, for untwisted and twisted fermionic fields, for several values of the ratio L/a. The 
dashed lines correspond to the fermionic condensate for the geometry of two plates in a space 
with topology R 4 , defined by Eq. f)3.25|) . Note that, for a massless field, the pure topological part 
vanishes. As it is seen from the graphs, for an untwisted (twisted) field the absolute value of the 
fermionic condensate increases with decreasing (increasing) length of the compact dimension. 



4 Vacuum expectation value of the energy- momentum tensor 

Another important local characteristic of the fermionic vacuum is the VEV of the energy- 
momentum tensor. In order to find this we use the mode-sum formula 

( T a-> = \ - (d^hu)^} , (4.1) 

P 

where the brackets denote the symmetrization over the indices enclosed. By taking into ac- 
count the expressions for the spinor eigenfunctions, we find the following expressions for the 
components of the vacuum energy-momentum tensor (no summation over /i) 

/r*\ = JMv f dkp V / M (An)/o; U2) 

2aV < nftzJ (2vr)^ 1 l-sin(2A n )/(2A n )' 1 ' 1 
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Figure 1: Fermionic condensate between two plates, in the model with spatial topology i? 3 x S , 
as a function of z/a for untwisted (04 = 0, left plot) and twisted («4 = 1/2, right plot) fields. 
The numbers near the curves correspond to the values of the ratio L^/a. The dashed lines in 
both plots correspond to the fermionic condensate for two plates on background of the space 
with trivial topology R A [see Eq. (|3.25j) ]. 

where u is defined in Eq. (|3.3|) and 

/ (0) (x) = w 2 [1 - cos (x) cos (:r(2z p+ 7a - 1))] , 

= -(kf/co 2 )f(°\x),l^0,p+l, (4.3) 
f {p+1) (x) = -kl +l . 

As in the case of the fermionic condensate, we assume the presence of the cutoff function in 
(14. 2ft . It can be checked that the VEVs of the separate components obey the trace relation 

{T£)=m{W), (4.4) 

with the fermionic condensate given by Eq. (|3.2p . 

The total vacuum energy (per unit volume along the directions z\, I = 1,2, ... ,p) in the 
region < z p+ \ < a, ^ Z[ ^ L[, I = p + 2, . . . , D, is obtained integrating (T§) over this region. 
By taking into account the expression for the energy density from Eq. (|4.2ft and Eq. (j2.9[) for 
the eigenvalues, we obtain the following expression for the vacuum energy 

n q eZi J V 7 n=l 

It coincides with the total energy of the fermionic vacuum in the region under consideration, 
evaluated as the sum of zero-point energies of elementary oscillators. Hence, we have shown 
that the vacuum energy obtained by integration of the energy density in the region between 
the plates does coincide with the energy evaluated as the sum of the zero-point energies of 
elementary oscillators. This means that the surface energy located on the boundaries is zero for 
the boundary conditions under consideration. Note that the surface energy vanishes for scalar 
fields with Dirichlet or Neumann boundary conditions as well, but this is not the case for a 
scalar field with Robin boundary conditions (see Refs. \29\ I32j). The total vacuum energy in 
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the region between the plates, defined by Eq. (|4.5|) . is investigated in Ref. |16j . Here we will be 
concerned with the vacuum energy density and stresses. 

For the evaluation of the VEVs, given by Eq. (|4,2|) . we apply to the sums over n the sum- 
mation formula (|3.4p . After steps similar to those already described in the case of the fermionic 
condensate, the VEV of the energy-momentum tensor can be expressed in the decomposed form 
(no summation over //) 



(Tp = (t^ 0) + (r;)« - -£-°5» 

Vq 

- f^-m^ )H)/2 



/ j / x+m 2ax i 1 



n g gZ<? 

where A p is given by Eq. (|3.16p and we have defined the functions 

x 2 ~ m L ( me 2xz me- 2xz 



W = — j— 1 ( 2 + — ,/i = 0,...,p, 

p + 1 V x — m x + m 1 



me 2xz me 2 ' 



x — m x + m 

2 



F^x) = k 2 [2 + '-^—^ - '——r~zr ) ' P = P + 2, . . . , -D, (4.7) 



iyx) = -2x 2 , /i = P + i. 

In Eq. (|4.6p . the term (no summation over ji) 



with /q ^ = k 2 +1 +m\^ /q = — /c 2 , / = 1, . . . ,p+ 1, is the VEV of the energy- momentum tensor 
for the topology R p+1 x (S 11 ) 9 when the boundaries are absent. This background is homogeneous 
and the corresponding densities are uniform. The second term on the rhs of Eq. (I4.6P is induced 
by the plate at z = when the second plate is absent (again no summation over jj): 



toV q J (2K) p n ^ z Jo JaP + yz + m^ W + ix m-ixj ' 

where 

fg»(z) = * 2 + k 2 + <, ft%)=0, 

fP(x) = —kf, 1 = 1,... ,p,p + 2,...,D, (4.10) 

The last term on the rhs of Eq. (|4.6p is induced by the presence of the second plate. Note that 
the vacuum stress along the direction normal to the plates vanishes in the geometry of a single 
plate and is uniform in the region between the two plates. The renormalized expressions for the 
purely topological part (T%) (°/ are given in Ref. |26| . In particular, for the corresponding energy 
density one has 

( T o ) (0) = 7jWE'^Vfi 9 ) 



f(D+i)/2(m^L 2 p+2 n 2 p+2 + ■■■ + L 2 D n 2 D ) 
: (L 2 p+2 n 2 p+2 + ---+L 2 D n 2 D )(^)/ 2 



(4.11) 
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where, as before, the prime means that the term n q = is excluded from the sum. Note that 
(no summation over fx) (Tjf)^ = (Tg)(°) for /j, = 1, . . . ,p + 1. In the discussion below we will 
be concerned with the boundary induced parts. 

The single plate part, Eq. (|4.9p . can be further simplified. With this aim, in the term with 
the exponent e 2lxz (e~ 2lxz ), we rotate the integration contour in the complex plane x by the 
angle tt/2 (— it/ 2). By using the integration formula (|3.14p . we find the following result (no 
summation over //) 



Cy)« _ - m MR V r ^-'^"Vm, (4-12) 

M K ^ J m (x + m)e 2xz v ; ' 



with the notations 

2 2 

F^(x) = — ^, /x = 0,l,...,p, 

p + 1 

^fr) = fi = p + 2,...,D, (4.13) 

and i ? Q P+1 ^(x) = 0. Note that the boundary induced parts in the vacuum stresses along the 
uncompactified directions parallel to the plates are equal to the boundary induced part in the 
energy density. This result is a consequence of the Lorentz invariance of the problem along these 
directions. The energy density is negative everywhere. In the absence of compact dimensions, 
from Eq. fj4. 12f) one finds (no summation over fi) 

(TpW=m(W)W/D, (4.14) 

for M = 0, ...D- 1, and (Tg)M = o. In Eq. (jOijl is given by Eq. (l3~T8|) . Of course, 

Eq. (|4.14p is a direct consequence of the trace relation (|4.4p . 

When the length of the j-th compact dimension is large, we replace the summation over 
rij in Eq. (|4.12p by an integral, with the help of Eq. (|3.19|) . As next step, we introduce a new 

integration variable v = ^jx 2 — . After changing to polar coordinates in the (y, v )-plane, 

the angular part of the integral is evaluated explicitly and, as it can be seen from Eq. (14. 121) . to 
leading order the result is obtained for the topology R p+2 x (S 1 ) q ~ 1 . 

For points near the plate, z <C L/, the dominant contribution to the series in Eq. (|4,12p 
comes from large values of ni. Replacing the corresponding summations by integrations, we see 
that, to leading order, the behavior of the VEV coincides with the one for topologically trivial 
space. If, in addition, z m~ 1 , combining Eq. (j4.14j) with Eq. (I3.20p . one finds (no summation 
over fi) 

, r mm- mN D T((D + l)/2) 

K {A^ D +^l 2 Dz D ' 1 ' 

for \i = 0, . . . j D — 1. For Lj<z the behavior of the VEVs is essentially different, depending on 
the phases ai in the periodicity conditions. When cti = 0, I = 0, . . . , D, the contribution of the 
nonzero modes is exponentially suppressed and the dominant contribution comes from the zero 
mode n q = 0. For the components along the uncompactified dimensions, to leading order we 
have (no summation over fj,) (T^)W w N D (Tji)^l +1 /(V q N p+1 ), fj, = 0, . . . ,p + 1, where 
is the corresponding VEV for two parallel plates in the space R p+1 located at z = and z = a. 
For the stress along the l-th compact dimension the contribution of the zero mode vanishes 
and, in the case under consideration (a; = 0), the dominant contribution comes from the modes 
with ri[ = ±1, rij = 0, j ^= I, and the stress is exponentially suppressed (no summation over I): 
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(T^)^ 1 ) oc e 2z / Ll . We have also assumed that L\ <C m. 1 . For 7^ the VEVs are exponentially 
suppressed and to leading order we have (no summation over fi) 

\ n> Wg(m z)(P+ 3 )/ 2 e 2m ° z ' 1 ' ' 

for /i = 0, . . . ,p and mo one recovers Eq. (|3.22p . We have a similar exponential suppression for 
the stresses along the compact dimensions. 

Now we return to the two plate geometry. Combining Eq. (|4.12p with Eq. (|4.6|) , for the VEV 
of the energy-momentum tensor in the region between the plates, we find (no summation over 

= to<°> -r&pY. r dx "^tr°'W- <4i7) 



where 



GJx) 



1 1 

P+l 



m 



2xz 1 2ax— 2xz\ 



x — m 



-(e AXZ + e 



M = 0,l,...,p, (4.18) 



G^{x) — k 2 



2 + ^—{e 2xz + e 1 
x — m 



V=p + 2,...,D, (4.19) 



and G p+ \{x) = F p+ \(x). From here, it follows that the boundary induced part in the vacuum 
energy density is always negative. It can be easily checked that the boundary induced parts 
obey the trace relation (|4.4p . For the VEV in the geometry of two parallel plates located at 
z = and z = a in a D-dimensional space with topology R D , one has (no summation over jj) 

Nn5 u r°° (r 2 - m 2 ^- 1 )/ 2 

w*> - - w^iW . * ss^ + i G * w - (420) 

where the expression for Gq^(x), with = 0, 1, . . . , D — 1, is obtained from Eq. (|4.18p by the 
replacement m nq — > m and Gqd(x) = —2x 2 . For a massless field this result reduces to 

Note that (T£) r d = (1 - 2- D )N D {T») { ^ , where (T^)^ is the corresponding VEV for a scalar 
field with Dirichlet or Neumann boundary conditions on the plates. 

For a massless field and in the presence of compact dimensions, the expressions for the VEV 
of the energy-momentum tensor are obtained from general formulas (|4.17p by putting m = 0. 
Alternative expressions are derived using the expansion (e y + 1) _1 = — X^rS=i( — l) n e~ ny ■ After 
integration one finds (no summation over /i) 

w> - to ,w + £ £ ° !?,(2o **» ) - (422) 

V / 9 n =l v ' n c ,eZi 



with the notations 



(x) = fp/2+1 0*0, M = 0, . . . , p, 

= (2an^) 2 / P / 2 (x), n = p + 2,...,D, (4.23) 



Ggi(x) = -(p + l)/ P / 2+ i(^)-x 2 / p/2 (^. 
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The corresponding vacuum densities are uniform. In this case the boundary induced part in the 
total energy (per unit volume along uncompactified dimensions) of the vacuum is aV q A(TQ) and 
the corresponding result obtained from Eq. (|4.22p coincides with the result derived in Ref. [16] 
by using zeta function techniques. 

Similar to Eq. (|3,27p . the VEV of the energy-momentum tensor may be presented in the 
decomposed form 

(TP = <r;}(°) + £ (T»)f + A(Tp, (4.24) 

3=1,2 

where {T£)*p is the part in the VEV induced by a single plate at z = aj, with a\ = 0, 02 = 0. 
The interference part in Eq. (I4.6P is given by the expression (no summation over /x) 



Vn 



„ ,°° (x 2 - m2 )(p-!)/ 2 



n„eZi " ln Q 



e 2ax + l 



(4.25) 



with the notation 



mt 



p+1 



m 



x + m 



m 



x + m 

-2xz 



e~ 2xz + e 



-2x(a—z) 



, fi = 0,1, 



+ e 



-2x(a— z) 



, V = P + 2 



, . . . , j 



D, 



,P, 



(4.26) 



and H p+ i(x) 



F p+ i(x). The surface divergences in the VEV of the energy-momentum tensor 



are contained in the single plate parts only and the interference part is everywhere finite. For 
a massless field, the single plate parts in the VEV of the energy- momentum tensor vanish and 
the interference part coincides with the second term on the rhs of Eq. (|4.22p . In the limit 
Li a and for a/ = the dominant contribution to the interference part comes from the zero 
mode rig. To leading order, we find (no summation over /j,) A(T^) = Ni)A{Tp RP +i / (y q N p+ i) , 
/i = 0, 1, . . . ,p+l, where A{Tp RP +i is the VEV for two plates in a (p+l)-dimensional spacetime 
with topology K p+1 . For the stress along the l-th compact dimension, the contribution of the 
zero mode vanishes. The dominant contribution comes from the modes with n\ = dbl, rij = 0, 



-2a/L t 



j 7^ I, and the stress is exponentially suppressed: A(T}) tx e 
assumed that Li -C m^ 1 . For a.[ 7^ and L\ <C a the interference part in the VEV of the 



where we have additionally 



energy-momentum tensor is suppressed by the factor e 



-2amo 



with mo defined in Eq. (I3.22I) . 



In the regions z p+1 < and z p+l > a the expression for the VEV of the energy-momentum 
tensor has the form 

TO = (t u J 0) + ro (1) > ( 4 - 27 ) 



where the boundary induced part is given by Eq. (|4.12p . with the replacements z — > \z\ and 
z — >■ z — a for the regions z < and z > a, respectively. In the presence of a constant gauge field 

, the formulas for the VEV of the energy- momentum tensor are obtained by the replacement 
(I2.12p . As in the case of the fermionic condensate, the vacuum energy- momentum tensor is a 
periodic function of A\L\ with the period of the flux quantum. 

The vacuum forces per unit surface of the plates are equal to the normal stress evaluated at 
z = and z = a. By taking into account that (T^^)^ = one finds, in the region between 
the plates, 



P(0+) = P(a-) = 
= -<2*S> (0) 



-<2*£>< > - AC/; 
2A P N D 



p+i 1 



E 

n„eZi 



ml )(f- x )/ 2 

11(7 / 



x—m 



+ 1 



(4.28) 
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where = X^£p+2p 7r ( n ' + ^z)/^/] 2 + m 2 - The term —{I^i}^ does not depend on the 
plate separation, it is a purely topological part in the vacuum forces. For this term, one has 
(Tp+l)^ = (T °)(°), with (T °}(°) given by Eq. EI) . The last term on the rhs of Eq. (PSD is 
induced by the presence of the second plate and determines the interaction forces between the 
plates. This term is negative and the interaction forces between the plates are always attractive, 
with independence of the periodicity conditions along compact dimensions and of the value of 
the gauge potential. In absence of the gauge field, Eq. (I4.28P coincides with the result obtained 
in Ref. [16], by differentiation of the total Casimir energy. For the vacuum forces in the regions 
z < and z > a, only the pure topological part contributes, and one has 

P(0-) = P(a+) = -(T;+ 1 )(°). (4.29) 

When the quantum field lives in all regions, the pure topological parts of the force acting 
from the left and from the right-hand sides of the plate compensate and the resulting force is 
determined by the last term on the rhs of Eq. (|4.28|) . In some important physical situations 
(bag model in QCD, finite-length carbon nanotubes, higher-dimensional models with orbifolded 
extra dimensions) the quantum field is confined to the interior of some region and there is no 
field outside. For the problem under consideration, if the quantum field is confined in the region 
between the plates, the total Casimir force, acting per unit surface of the plate, is determined 
by Eq. (|4.28|) and the pure topological part contributes as well. The resulting force can be either 
attractive or repulsive, depending on the phases in the periodicity conditions along the compact 
dimensions, and also on the gauge potential. In this case, the Casimir effect could be used as a 
stabilization mechanism for both the interplate distance and the size of the compact subspace 
in Kaluza-Klein-type models and in braneworlds. This is a quite remarkable and very useful 
result, because this force could be, in principle, easily controlled. 

In Fig. we have plotted the ratio of the boundary induced part in the Casimir energy 
(per unit surface along uncompactified dimensions) for two parallel plates in the spacetime with 
topology R 3 x S 1 by the Casimir energy in R 3 , E$ xS1 /E R3 = L{T$)W /(T$} R3 (L being the 
length of the compact dimension), with (Tq)^ = (Tq) — (Tq)(°\ for a massless fermionic field, 
as a function of a/L. The values on each of the curves correspond to those of the parameter 
q 4 . Note that from Eq. (OT1) one has (T °) R 3 = -7vr 2 /(2880a 4 ). The feature described before 
is apparent: for large values of a/L the Casimir energy is suppressed for ^ 0. 

For a massless case the fermionic field is conformally invariant. We can generate the cor- 
responding VEVs in conformally flat spacetimes by using standard conformal transformation 
techniques (see, for instance, |28l ) and the formulas given above. Consider a conformally-flat 
spacetime with the line-element 

D 

ds 2 = n 2 (zi)(dt 2 -^2(dz l ) 2 ), (4.30) 
t=l 

and spatial topology R p+1 x (iS 1 ) 9 , with ^ zi ^ L\, I = p + 2, . . . , D. As before, we assume the 
presence of two boundaries located at z = and z = a = z), with the boundary conditions 

[1 + ^7(Q) n (n)/i]^(n) = 0> where the subscript $7 specifies the quantities on the background 
described by Eq. (|4.30p . For curved space gamma matrices we have 7^ = ef^ 1 , with the tetrad 
field = Q~ 1 5j l . By taking into account that, under the conformal transformation for a massless 
field, one has = ^ _D V' an d n (Q)fi = Qn^, we see that the MIT bag boundary condition 
is conformally invariant. For the fermionic condensate and the VEV of the energy- momentum 
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2.0 




Figure 2: Ratio of the boundary induced part in the Casimir energy for two parallel plates in the 
spacetime with spatial topology R 3 x S 1 by the standard Casimir energy in i? 3 , for a massless 
fermionic field, as a function of a/L. The values on each of the curves correspond to those of 
the parameter a^. 

tensor on the spacetime (I4.30p . one finds 

TO(o) = (Tp { n), R o+n- D - l (Tp, (4.31) 

where (^V')(n),i? 13 an d are the corresponding VEVs for spacetime (I4.30P with trivial 

spatial topology R D . Note that for points away from the boundaries, all divergences are con- 
tained in these terms and the renormalization procedure is needed for them only. The second 
terms on the rhs of Eq. (|4.3ip are induced by the non-trivial topology and by the boundaries. 
Most important special cases of Eq. (|4.30p are the de Sitter (dS) and the anti-de Sitter (AdS) 
spacetimes, with f2jL = l/(Ht) 2 and fi^dS = l/(kz) 2 , described in the inflationary and Poincare 
coordinates, respectively. In particular, the results for AdS bulk can be applied to massless Dirac 
fields in higher-dimensional Randall-Sundrum-type braneworld models with two branes and with 
compact internal spaces. 

5 Casimir densities in carbon nanotubes 

In a significant variety of planar condensed matter systems, the low-energy sector is very well 
described by the Dirac-like model. A well-known example is the important case of graphene. In 
this section we apply general results obtained above for electrons in cylindrical nanotubes of finite 
length. Recently, carbon nanotubes have attracted a lot of attention due to the experimental 
observation in them of a number of novel electronic properties, what renders them very important 
for technological purposes. A single-wall cylindrical nanotube is a graphene sheet rolled into a 
cylindrical shape. The low-energy excitations of the electronic subsystem in a graphene sheet 
can be described by a pair of two-component spinors, tpA and ipB, corresponding to the two 
different triangular sublattices of the honeycomb lattice of graphene (see, for instance, [2J[3]). 
The Dirac equation for these spinors has the form 

(iVp^Do + vy l Di - m)ipj = 0, (5.1) 



17 



where J = A, B, I = 1, 2, and = <9 M + ieA^ with e = — |e| for electrons. In Eq. (|5.1|) . 

~ 10 8 cm/s represents the Fermi velocity which plays the role of the speed of light. To make 
the treatment more general, we have included in Eq. (|5.ip the mass (gap) term. The gap in the 
energy spectrum is essential for many physical applications. It can be generated by a number 
of mechanisms (see, for example, [3 [331 EH S3 [36] ) . I n particular, they include the breaking of 
symmetry between two sublattices by introducing a staggered onsite energy [3] , the phenomenon 
of magnetic catalysis [33J , and the deformations of bonds in the graphene lattice |34| . Another 
approach is to attach a graphene monolayer to a substrate, the interaction with which breaks the 
sublattice symmetry [35]. Note that the Casimir interaction between graphene sheets, resulting 
from the quantum fluctuations of the bulk electromagnetic field, has been recently investigated 
in Ref. [37] (for the comparison of the results based on the hydrodynamic and Dirac models of 
dispersion for graphene, see Ref. [38]). 

For the geometry of a carbon nanotube the spatial topology is R 1 x S . The nanotube is 
characterized by its chiral vector Ch = n w ai + m w a2, where n w , m w are integers, ai and &2 are 
the basis vectors of the hexagonal lattice of graphene and a = |ai| = |&2 1 = 2AQA is the lattice 
constant. For the length of the compact dimension, one has L = \Ch\ = t}yn| + + n w m w , 
where for zigzag and armchair nanotubes, Ch = (n w ,0) and Ch = (n w ,n w ), respectively. All 
other cases correspond to chiral nanotubes. In the case n w — m w = 3q w , q w G Z, the nanotube 
will be metallic and in the 7^ 3q w the nanotube will be a semiconductor with an 

energy gap inversely proportional to its diameter. We will assume that the nanotube has finite 
length, a. As the Dirac field lives on the cylinder surface, it is natural to impose bag boundary 
conditions (12. 2p on the cylinder edges which insure a zero fermion flux through these edges. 
The additional confinement of the fermionic field along the tube axis leads to the change of the 
VEVs. The corresponding expressions for the fermionic condensate and the energy-momentum 
tensor are obtained from the formulas of the previous sections, taking D = 2, p = 0, q = 1. The 
periodicity condition along the compact dimension for the fields ipj depends on the chirality of 
the nanotube. For metallic nanotubes, we have periodic boundary condition (a; = 0) and for 
semiconductor nanotubes, depending on the chiral vector, there are two classes of inequivalent 
boundary conditions, corresponding to a\ = ±1/3. These phases have opposite signs for the 
sublattices A and B. The presence of the gauge field in Eq. (|5.ip leads to an Aharonov-Bohm 
effect in carbon nanotubes [39]. This effect manifests itself in a periodic energy gap modulation 
and conductance oscillations, as a function of the enclosed magnetic flux, with a period of the 
order of the flux quantum. As we will see below, similar oscillations arise for the fermionic 
condensate and Casimir densities. 

First, we consider the fermionic condensate (see Ref. [40J for the fermionic condensate and 
VEV of the fermionic current in a (2 + l)-dimensional conical spacetime, in the presence of 
a circular boundary). The corresponding expression is obtained from the formulas in Sect. [3] 
taking p = 0, D = 2 and summing the contributions coming from the two sublattices with 
opposite signs of ai = a. For an infinite carbon nanotube, from Eq. (|3. 10[) one has, for the pure 
topological part, 

2mv F — <>- nmL 



M)® = V cos(2vrna) cos(27rn$/$ ), (5.2) 

ttL ' n 

n=l 

where $ = A2L is the magnetic flux across the nanotube and a = 0, 1/3 for a metallic and for a 
semiconducting nanotube, respectively. In the case of a nanotube of finite length, a, the general 
expression (|3.23p takes the form 

e r^;""-^-^, (5 .3) 

n=-oc ]=+,- Jm "' ^j x 2 _ m 0) 2 ( f±^ e 2a* + 1 
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with the notation 

= \J k { n )2 + m 2 , fcW = 2Tr\n + a± $/§ \/L. (5.4) 

Note that for metallic nanotubes the contribution of the terms with j = + and j = — coincide. 
In the absence of the magnetic flux, the sublattices give the same contribution to the condensate. 
For a massless field the purely topological part vanishes and the boundary induced part reduces 
to 

_. +00 00 

W)- = ^Z E E W £(-l)'[*i(2*W(ai - *)) + K^ial - a + z))]. (5.5) 

n=—ooj=+,— 1=1 

As already mentioned, the condensate is a periodic function of the magnetic flux, with period 
equal to the flux quantum <3?o- Various important limiting cases directly follow from the analysis 
given above. 

Now we turn to the VEV of the energy-momentum tensor. For the pure topological part, 
we have the expression (no summation over /j,) 

2 VF § U 00 e -nmL 

( T DcJ = —TT 1 cos(2vrna) cos(2vrn$/$ )C. (nmL) — , (5.6) 

n=i 

with the notations 

C (x) = d(x) = l + x, C 2 (x) = -2-2x- x 2 . (5.7) 

In particular, in absence of magnetic flux, the corresponding energy density is positive for 
metallic nanotubes and negative for semiconducting ones. This means that, from the topolog- 
ical viewpoint, semiconducting nanotubes are more stable. For a finite length nanotube, from 
Eq. (|4.17p we find the following expression (no summation over jj) 



VT? W +°° roc (Jl (j')2x_i/2 

TO- - - -§ E E mU) {x ^:L ] +1 oyw. (") 

n=-oo j=+- m ™ x—m 



where 



G^\x) = (x 



n 



2 -|- — i y c 2xz -\- £>2ax—2xz\ 

x — m 



G ( i\x) = -2x 2 , (5.9) 



2 -|- — (c 2xz -\- e 2ax ~ 2xz ] 



m 



In absence of magnetic flux, the energy density corresponding to Eq. (|5.8|) is always negative 
for semiconducting nanotubes. For metallic nanotubes the energy density is positive for long 
tubes and negative for short ones. The forces acting on the tube edges are determined by the 
component (T^) cn . As already explained before, the force corresponding to the second term on 
the rhs of Eq. (j5.8[) is attractive, with independence of the nanotube chirality and of the magnetic 
flux. This term dominates for short nanotubes. As regards the first term, it dominates for long 
tubes and, in the absence of magnetic flux, it corresponds to the attractive force for metallic 
nanotubes and to the repulsive force for semiconducting ones. Hence, in absence of a magnetic 
flux, the resulting force is always attractive for metallic nanotubes, whereas for semiconducting 
ones the force is attractive for short tubes and repulsive for long ones. The sign of the force for 
long nanotubes can be simply controlled by tuning up the magnetic flux. This result is again of 
high technological importance. 
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In the massless case, the above formulas take the form (no summation over fi) 

TO™ = TOg + E 4 0) (2-4 j) )> (5-10) 

s=l j'=+,— n=— oo 

with the notations 

G[, 0) (x) = xKxix), G { ®\x) = x 2 K (x), 

Gf\x) = -xK l {x)-x 2 K (x). (5.11) 

In Fig. [3] we have plotted the vacuum energy density and vacuum stresses for a massless fermionic 
field in metallic nanotubes as functions of the tube length in absence of magnetic flux, ao being 
an arbitrary length scale. 




a/a a/a 



Figure 3: VEV of the energy density (left plot) and vacuum stresses (right plot) for a metallic 
nanotube, as functions of the tube length, for several values of its circumference, L/ao (numbers 
near the curves). On the left plot full/dashed curves correspond to / = l/l = 2, respectively. 

For long tubes the pure topological part dominates in the VEV of the energy density and in 
the case of the values for parameters corresponding to Fig. [3] one has: 

(Tl)cn « (T/)W = ^Mdiag(l, 1, -2), (5.12) 

with £(3) ~ 1.202. The effective pressure on the tube edges is defined by —(T^) cn and, as seen 
from the right plot, it is always negative, corresponding to attractive forces between the edges 
for metallic nanotubes. 

Fig. H] is a corresponding plot of the VEVs of the energy density and stresses for semicon- 
ducting nanotubes (a = 1/3), for a massless field in the absence of magnetic flux. Now, for long 
tubes we have the limiting value 

(T}U - (Tf)S = -^^diag(l, 1, -2). (5.13) 

The forces acting on the tube edges are attractive at small distances and repulsive at large 
distances. 

In Fig. [5] we have plotted the vacuum energy density as a function of the magnetic flux 
for metallic (left plot) and semiconducting (right plot) nanotubes, respectively. Recall that 
the vacuum densities are periodic functions of the magnetic flux, with period equal to the flux 
quantum. The numbers labelling the curves correspond to the values of the ratio L/a. 
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Figure 5: VEV of the energy density for metallic (left plot) and semiconducting (right plot) 
nanotubes, as functions of the magnetic flux in units of the flux quantum. The curves correspond 
to different values of the ratio L/a (labels of the curves). 
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6 Conclusions 



In this paper we have considered the effect of compact spatial dimensions on the fermionic con- 
densate and VEV of the energy-momentum tensor for a massive fermion field in the geometry 
of two parallel plates on which the field obeys MIT bag boundary condition. Along the com- 
pact dimensions, we have assumed periodicity conditions (|2.ip with constant phases a/. The 
eigenvalues of the wave-vector component normal to the plates are roots of the transcendental 
equation (12. 9p . The mode sums for the fermionic condensate and the energy- momentum ten- 
sor contain series over these eigenvalues. By applying analytic continuation techniques, as the 
Abel-Plana-type summation formula, to these series, we have been able to explicitly extract and 
separate, in a cutoff independent way, the purely topological part and the contributions induced 
by the single plates. Purely topological contributions were investigated in Ref. [26], and here we 
have been mainly concerned with the boundary induced parts. 

In a (-D + l)-dimensional spacetime with spatial topology R p+1 x (S 1 ) q and with two parallel 
boundaries, the fermionic condensate is given by Eq. (13.230 . there the first term on the rhs is 
a purely topological contribution and the second term is induced by the presence of the plates. 
The boundary induced part is always negative, whereas the purely topological one, given by 
Eq. (|3.10p . can be either positive or negative, depending on the values of the phases in the 
periodicity conditions along the compact dimensions. For a massless field, the general formula 
for the fermionic condensate could be further simplified to Eq. (|3.24p . Extracting there the 
contributions corresponding to the geometry of a single plate, the fermionic condensate could 
also be presented in the alternative form (|3.27p . The fermionic condensate induced by a single 
plate, located at z = 0, is given by Eq. (|3.17p . The condensate diverges on the boundary. 
For points near the boundary, the leading term in the asymptotic expansion over the distance 
from the plate is given by Eq. (|3.20p . This term does not depend on the lengths of the compact 
dimensions and coincides with the boundary induced contribution of the fermionic condensate for 
a single plate, in a space with trivial topology R D . Far from the plate, the asymptotic behavior of 
the fermionic condensate essentially depends on the phases present in the periodicity conditions 
along the compact dimensions. 

For ai = 0, I = p + 2, . . . , D, the main contribution comes from the zero mode and the 
quantity V q (^ip)p^/N]j coincides with the corresponding result for a plate in topologically trivial 
(p+ l)-dimensional space, R p+1 . For ai ^ 0, the zero mode is absent and the boundary induced 
part in the fermionic condensate is suppressed by the factor e ~ 2m ° z , where mo is defined in 
Eq. (|3.22p . The interference part in the representation (13.271) for the fermionic condensate 
is given by Eqs. (|3.28[) and f|3.29[) for massive and for massless fields, respectively. Surface 
divergences are contained in the single plate parts and the interference contribution is finite and 
positive everywhere. 

The VEV of the energy-momentum tensor was investigated in Sec. HI In the region between 
the plates, this VEV is given by Eq. (|4.17j) . where the second term on the rhs is induced by the 
plates. The vacuum stresses along the uncompactified dimensions parallel to the plates are equal 
to the energy density. The boundary induced contribution to the energy density and stresses 
along the compact dimensions are always negative, whereas the corresponding contribution to the 
stress normal to the plates is positive and, moreover, a uniform function. For a massive field, the 
energy density and stresses along the directions parallel to the plates depend on the distance from 
the plate and diverge on the boundaries. For a massless field the VEV of the energy-momentum 
tensor is uniform, and given by Eq. (|4.22p . An alternative representation for the VEV of the 
energy- momentum tensor is (I4.24|) . with the interference part defined in Eq. (|4.25|) . The single 
plate contributions, in this decomposition, are given by (|4.12p . The surface divergences in the 
VEVs are contained in the single plate parts and the interference part is finite everywhere. Near 
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the plates, the VEV of the energy-momentum tensor is dominated by the single plate parts 
and the leading term in the corresponding asymptotic expansion is given by Eq. (|4.15p . In the 
geometry of a single plate and at large distances from the plate, compared to the length of the 
compact dimension, the behavior of the VEV is essentially different and crucially depends on 
the phases a\ included in the periodicity conditions. Thus, for periodic boundary conditions, 
ai = 0, the dominant contribution comes from the zero mode and for the components along the 
uncompactified dimensions one has (no summation over fi) (Ttf)^ Nd{T^)^ +1 / '(V q N p+ i), 
/i = 0, . . . ,p + 1, where (T^)^ p+1 is the corresponding VEV for plates in the space R p+l . The 
stress along the l-th compact dimension is suppressed by the factor e~ 2z l Ll . For ct\ ^ all 
components of the vacuum energy-momentum tensor are suppressed by the factor e ~ 2rn o z . 

The formulas derived for the fermionic condensate and for the VEV of the energy-momentum 
tensor are easily generalized to the case when a constant gauge field is present. In spite of the 
fact that the corresponding magnetic field strength vanishes, the non-trivial topology of the 
configuration leads to an Aharonov-Bohm-type effect on the vacuum expectation values. The 
corresponding formulas are obtained from those in the absence of a gauge field by the replacement 

(1232]). 

The vacuum force per unit surface of the plates is determined by the normal stress and, 
for the region between the plates, it is given by Eq. (|4.28p . where the first term is the pure 
topological part and the second one the interaction part of the force. When the quantum field 
lives in all regions of space, the pure topological parts of the force acting from the left- and 
from the right-hand sides of every plate compensate each other, and the resulting force is just 
determined by the interaction part. This turns out to be always attractive, with independence 
of the periodicity conditions along the compact dimensions and the value of the gauge potential. 
In more physical situations, however, when the quantum field is confined in the region between 
the plates, the pure topological part contributes to the resulting force as well and in these cases 
the final forces can be either attractive or repulsive. Remarkably, the Casimir effect can then 
be used as a stabilization mechanism for both the interplate distance and for the size of the 
compact subspace in Kaluza-Klein-type models and in braneworld theories. 

For massless fermionic fields, the Casimir densities for the two-plate geometry in conformally- 
flat spacetimes with compact spatial dimensions are obtained from the results above simply by 
using standard conformal transformation techniques. The corresponding fermionic condensate 
and the VEV of the energy-momentum tensor are thus given by Eq. (|4.3ip . where the separate 
contributions: the one induced by non-trivial topology and the other by the boundary condi- 
tions, are given by the second terms on the right-hand sides. The most important examples of 
configurations of this type are dS and AdS spacetimes. 

In Sect. [5j we have applied our general results to the electrons in finite- length graphene 
nanotubes. The long-wavelength excitations of this system is basically described by Dirac's 
theory, with the Fermi velocity playing the role of speed of light. The corresponding expressions 
for the fermionic condensate and for the Casimir densities are obtained from the results of 
the preceding sections by summing the contributions of the two triangular sublattices of the 
honeycomb lattice of the graphene sheet, with opposite signs of the phases along the compact 
dimension. In the presence of a magnetic flux, the fermion condensate is given by Eqs. (|5.3p 
and (15. 5p . for massive and for massless fields, respectively. The pure topological contribution is 
then given by (|5.2p . In these formulas, the values a = and a = 1/3 correspond to metallic and 
to semiconducting nanotubes, respectively. The corresponding expressions for the VEV of the 
energy-momentum tensor have the form (|5.8p and (|5.10p . The VEVs are periodic functions of 
the magnetic flux with period equal to the flux quantum. In the absence of magnetic flux, the 
Casimir forces acting on the edges of the nanotube are always attractive for metallic nanotubes. 



23 



However, for semiconductor nanotubes the forces are attractive for short tubes and repulsive for 
longer ones. In the presence of the magnetic flux, the sign of the force for long tubes can be 
controlled by tuning the flux. This possibility opens the way to the design of efficient actuators 
driven by the Casimir force at the nanoscale, what is a most cherished goal of present day 
technology. 
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